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NANOSTRUCTURED  COMPOSITES:  EFFECTIVE  MECHANICAL  PROPERTY 
DETERMINATION  OF  NANOTUBE  BUNDLES 

E.  Saether,*  R.  B.  Pipes1,  S.J.V.  Frankland11’** 

Abstract.  Carbon  nanotubes  naturally  tend  to  form  crystals  in  the  form  of  hexagonally  packed 
bundles  or  ropes  that  should  exhibit  a  transversely  isotropic  constitutive  behavior.  Although  the  intratube 
axial  stiffness  is  on  the  order  of  1  TPa  due  to  a  strong  network  of  delocalized  bonds,  the  intertube  which 
is  controlled  by  weak,  nonbonding  van  der  Waals  interactions  is  orders  of  magnitude  less.  An  accurate 
determination  of  the  effective  mechanical  properties  of  nanotube  bundles  is  important  in  order  to  assess 
potential  structural  applications  such  as  reinforcement  in  future  composite  material  systems.  A  direct 
method  for  calculating  effective  material  constants  is  developed  in  the  present  study.  The  Lennard-Jones 
potential  is  used  to  model  the  nonbonding  cohesive  forces.  A  complete  set  of  transverse  moduli  are 
obtained  and  compared  with  existing  data. 

Key  Words,  carbon  nanotubes,  elastic  constants,  molecular  simulation 
Subject  Classification.  Materials  Science 

1.  Introduction.  Future  nanostructured  composite  materials  are  expected  to  incorporate  carbon 
nanotube  reinforcement  either  dispersed  individually  or  as  nanofilamentary  bundles  or  ropes  yielding 
unprecedented  mechanical  properties.  A  carbon  nanotube  is  a  cylindrical  molecule  composed  of  single  or 
multiple  walls  of  graphene  sheets.  These  sheets  are,  in  turn,  composed  of  hexagonal  units  or  graphene 
rings  of  carbon  atoms  that  are  bonded  through  highly  stable  sp~  hybridized  orbitals.  A  typical  carbon 
nanotube  is  schematically  depicted  in  Figure  1  while  Figure  2  contains  a  cross-section  of  a  bundle 
ensemble  of  individual  nanotubes  obtained  through  transmission  electron  microscopy  (TEM)[1], 

Numerous  analytical  and  experimental  studies  have  been  made  to  determine  the  elastic  properties  of 
individual  nanotubes.  Axial  Young’s  moduli  on  the  order  of  1  TPa  have  been  measured  using  atomic 
force  microscopy  (AFM)[2]  and  thermal  vibrations  [3].  Analytical  studies  have  utilized  ab  initio 
calculations  [4],  tight-binding  methods  [5],  molecular  dynamic  simulations  (MD)[6]  and  lattice 
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dynamics[7,8].  These  elastic  properties  are  entirely  based  on  the  strong  intratube  valence  forces  of  the 
carbon-carbon  bonds. 


Cl„,mr ,  c-  i  ii  j  FIGURE  2.  Typical  nanotube  bundle [1], 

FIGURE  1 .  Single-walled  nanotube. 


Nanotube  ensembles,  however,  typically  form  hexagonally  packed  crystal  configurations  in  which  the 
intertube  force  interactions  are  due  primarily  to  non-bonding  van  der  Waals  effects  which  are  much 
weaker  than  the  valence  forces  and  are  highly  nonlinear.  Less  consideration  has  been  given  to  the 
transverse  mechanical  properties  of  nanotube  bundles  which  depend  on  a  good  description  of  these  non- 
bonded  interactions.  These  intertube  cohesive  properties  are  of  special  interest  for  use  in  predicting  the 
properties  of  carbon  nanotube  polymer  composites[9]  and  fibers  of  woven  nanotubes [10].  Selected 
moduli  of  nanotube  ropes  have  been  calculated  with  a  continuum  model  based  on  the  integrated  average 
of  the  discrete  Lennard-Jones  potential[ll],  MD  simulation  using  the  Tersoff-Brenner  potential[12]  and 
lattice  dynamic  methods [13, 14], 

For  the  present  study,  a  direct  summation  of  atom-pair  potentials  is  used  to  avoid  any  simplifications 
made  to  the  nonlinear  van  der  Waals  interactions  and  to  permit  optimum  flexibility  in  representing 
discontinuous  or  curved  nanotubes.  Additionally,  a  direct  method  permits  irregular  and  disrupted 
nanotube  lattices  to  be  modeled.  Because  the  fundamental  constituents  of  nanotube  bundles  are  only 
resolvable  at  nanometer  length  scales,  analyses  to  predict  macroscopic  properties  must  necessarily  merge 
concepts  and  techniques  from  continuum  elasticity  theory  and  discrete  molecular  simulation.  The  basic 
approach  of  subjecting  a  molecular  ensemble  to  applied  strain  modes  and  recovering  effective  moduli 
from  energy  measures  has  been  used  in  molecular  dynamic  simulations [15, 16].  The  methodology 
developed  herein  combines  a  unit  cell  continuum  model  with  molecular  static  calculations  to  determine 
effective  moduli  in  aligned  carbon  nanotube  bundles. 
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The  Lennard-Jones  potential  is  utilized  to  simulate  the  van  der  Waals  interaction  forces  among  carbon 
atom-pairs  in  aligned  carbon  nanotube  arrays.  An  achiral  “zig-zag”  configuration  is  assumed  for  the 
carbon  nanotubes  and  the  tube  radius  is  assumed  small  such  that  the  cross-section  can  be  considered  rigid. 

The  objective  of  this  work  is  to  formulate  a  molecular  mechanics  model  to  understand  and  predict  the 
elastic  properties  of  carbon  nanotube  bundles.  This  method  utilizes  a  unit  cell  approach  for  determining 
selected  transverse  moduli  of  a  hexagonal  crystal  of  aligned  nanotubes.  The  predicted  moduli  are  shown 
to  exhibit  the  transverse  isotropy  anticipated  for  a  material  possessing  hexagonal  symmetry.  The  predicted 
moduli  are  compared  with  available  published  data. 

2.  Modified  Unit  Cell  Formulation.  In  elasticity  analyses,  the  method  of  unit  cells  has  been 
developed  to  determine  the  effective  properties  of  heterogeneous  materials  by  identifying  and  analyzing 
convenient  domains  of  repeating  microstructure.  Applied  to  the  determination  of  effective  continuum 
elastic  moduli  of  bundles  of  aligned  carbon  nanotubes,  a  minimal  repeating  unit  of  nanotubes  is  defined 
and  subjected  to  continuous  field  deformation  modes.  An  assemblage  of  a  primary  unit  cell  of  nanotubes 
with  surrounding  image  regions  that  are  required  in  applying  periodic  boundary  conditions  (pbc)  is 
presented  in  Figure  3.  In  this  approach  the  nanotubes  are  regarded  as  sources  of  potential  energy  that 
require  interactions  with  nanotubes  outside  the  unit  cell.  Under  periodic  boundary  conditions,  cells  of 
nanotubes  in  the  transverse  plane  and  nanotube  segments  in  the  axial  dimension  are  treated  as  images  of 
the  constituents  within  the  cell  and  used  in  the  calculation  of  potential  energy.  This  permits  interactions 
between  atom-pairs  across  the  boundary  to  avoid  introducing  discontinuities  in  the  force  field.  In  general, 
these  conditions  ensure  conservation  of  mass  and  energy,  avoid  surface  or  boundary  effects,  and 
mathematically  give  the  primary  unit  cell  a  strict  periodicity  such  that  it  can  be  considered  to  represent  an 
infinite  ensemble  of  molecules. 

By  combining  concepts  from  continuum  elasticity  and  molecular  dynamics,  these  representative  units 
will  be  referred  to  herein  as  ‘pbc-unit  cells’.  Pbc -unit  cells  can  be  constructed  of  arbitrary  order  but,  with 
a  proper  definition  of  the  repeat  geometry,  the  unit  energy  of  the  primitive  cell  remains  the  same. 
Therefore  the  lowest  order  cell  is  used  for  computations.  A  minimum-order  hexagonal  pbc-unit  cell  is 
shown  in  Figure  4. 

The  initial  equilibrium  configuration  of  the  hexagonal  unit  cell  is  determined  by  minimizing  the 
energy  of  the  system  as  the  nanotubes  are  moved  radially  outward  from  a  fixed  center.  This  establishes 
the  equilibrium  radius,  Req,  and  the  nanotube  center-to-center  separation  distance,  S,  as  shown  in  Figure  5. 

A  rigorous  definition  of  the  pbc-unit  cell  dimensions  is  required  to  ensure  invariance  of  the  unit 
energy  with  cell  size.  The  required  planar  area  of  the  pbc-unit  cell  is  given  by 


3 


(2.1) 


ACE1I  =|m/3r; 

where  M  is  the  number  of  nanotubes  within  the  unit  cell.  The  effective  depth  of  the  pbc-unit  cell  is 
obtained  by  first  selecting  a  number  of  repeat  units  (circumferential  rings  of  graphene).  Kseg,  and  adding 
one  additional  unit  to  account  for  image  segments  in  the  axial  dimension.  The  total  area  of  the  enclosed 
graphene  units  is  then  equated  to  the  surface  area  of  a  perfect  cylinder  having  the  same  radius  as  the 
nanotube  which  is  expressed  as 

(2-2)  N  Agu  (Kseg  + 1)=  2  JtRnl  deff 

where  N  is  the  number  of  graphene  units  around  the  nanotube  circumference,  Agu  is  the  area  of  an 


b)  (l,2)-plane  section 

FIGURE  3.  Pbc-unit  cell  showing  outside  periodic  image  nanotubes  within  the  (2,3)- 
plane 


FIGURE  4.  Minimal  hexagonal  pbc-unit  cell 
with  surrounding  image  cells. 


FIGURE  5.  Equilibrium  radius  definition 
for  a  hexagonal  cell. 


individual  graphene  unit,  and  R„r  is  the  nanotube  radius.  The  radius  of  an  achiral  nanotube  can  be  derived 
as 
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(2.3) 


V3 

Rnt=YbSec 


L2v 


1-  — 
N 


where  b  is  the  carbon-carbon  bond  length.  With  the  area  of  a  graphene  unit  given  by  Agu  =  3V3b2/2, 
Equation  (2.2)  yields  the  effective  depth  of  the  pbc-unit  cell  as 


(2.4) 


dt£f  =  — NbCos 
n 


L-2  v 


1-  — 

N 


(k,8+i) 


3.  Material  Constitutive  Relationship.  The  force  field  within  a  nanotube  crystal  consists  of  a 
combination  of  strong  linear  bonding  forces  acting  within  the  nanotube  and  weak  non-bonding  forces 
acting  between  adjacent  nanotubes.  This  disparity  between  the  magnitude  of  interatomic  forces  leads  to  a 
highly  anisotropic  constitutive  relation.  The  hexagonal  symmetry  of  the  nanotube  crystal  dictates  a 
transversely  isotropic  stress-strain  relation  given  by 

G  =C  8  +C  8  +  C  8 

W  11  v"ll  °11  1  '-'12  °22  1  '-'12  °33 

G  =C  8  +  C  8  +  C  8 

W  22  '-'12  °11  ^  '-'22  °22  ^  ^23  °33 


(3.1) 


G  =  C  8  +  C  8  +  C  8 

W  33  V"'  12  °11  ^  ^  23  °  22  ^  '-'  22  °  33 


°.2  =2C„e12 


=2C44e13 


*^23  (G22  ^23)^23 


Additional  relationships  between  the  material  constants  are  given  in  Reference  17  for  a  hexagonal  system 
as 


(3.2) 


K23=  (C22+C23)/2 
G23  =  (C22  -C23)/2 


This  material  is  uniquely  defined  by  five  independent  quantities  where  the  ‘  1  ’  direction  is  directed  along 
the  tube  axis.  The  current  effort  will  focus  on  the  prediction  of  normal  and  shear  properties  in  the 
transverse  plane. 
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4.  Analysis  Methodology.  The  methodology  used  to  determine  selected  nanotube  crystal  properties 
involves  defining  an  appropriate  pbc-unit  cell,  applying  selected  strain  modes  to  the  crystal,  and 
computing  the  potential  energy  due  to  atom-pair  interactions  as  a  function  of  the  deformation  kinematics. 
Effective  elastic  constants  are  then  determined  from  the  variation  in  the  system  potential  energy  as 


(4.1) 


c,= 


vv. 


C)  2  <t> 

3ei  3e. 


where  Qj  is  the  material  stiffness,  V  is  the  pbc-unit  cell  volume,  <f>  is  the  potential  energy,  and  £k  is  an 
applied  strain  mode. 

Strain  modes  are  applied  to  the  nanotubes  in  the  crystal  by  the  imposition  of  specific  deformation 
fields.  The  G23  shear  modulus  for  a  hexagonally  packed  nanotube  array  is  calculated  using  a  pbc-unit  cell 
subjected  to  a  pure  shear  strain  mode  as  shown  in  Figure  6.  The  magnitude  of  the  shear  strain  is  given  by 
twice  the  shear  angle  or  y23  =  2  0  .  During  a  progressive  deformation  with  increasing  0,  the  potential 
energy  is  computed  by  summing  all  atom-pair  interactions  between  adjacent  nanotubes.  The  G23  shear 
modulus  is  then  obtained  from  the  elastic  strain  energy,  <t> ,  using  a  finite  difference  approximation  as 

(4.2)  G  =  1  30  = _ Az° _ 

23  Vol3Y;3  A^cUAjJ2 

where  Acel]  is  the  planar  area  of  the  pbc-unit  cell,  dee  is  the  effective  length  of  the  nanotubes,  Ai  and  A2  are 
first  and  second  order  central  difference  operators,  respectively  . 
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FIGURE  6.  Imposed  shear  deformation  on  hexagonally  packed  nanotube  array. 
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The  bulk  modulus  is  computed  by  applying  a  dilatational  strain  as  shown  in  Figure  7.  Because  the 
strain  in  the  axial  dimension,  £n,  is  assumed  to  be  zero,  the  dilation  is  defined  as  e  =  £22+£33  with 
e22  =  £ 33  =  £ .  The  bulk  modulus  is  obtained  by  applying  Equation  (7)  using  a  strain  given  by  2e.  The 
calculation  of  the  E22  Young’s  modulus  and  the  Poisson’s  ratio  v 23  is  performed  by  applying  £22  strain 
increments  in  the  2-direction  and  repositioning  the  tubes  in  the  3-direction  to  minimize  the  energy.  The 
transverse  repositioning  of  the  tubes  perpendicular  to  the  load  axis  directly  gives  a  measure  of  the  £23 
strain  from  which  the  Poisson  ratio  is  determined  as  n23  =  -e23  /e22  .  This  is  depicted  in  Figure  8. 


x3 

A 


Q.Q.Q.O 

OOOOO 

OOOOOO 

OOOOOOO 

OOOOOO 

00  OQO 
0000 


>  X2 


0000 
o/o o"o\o 
0/0  o  o  o\o 
o<o  000  o>o 
o\o  o  o  0/0 
o\o  o  0/0 
o'o'o'o 


FIGURE  7.  Imposed  dilatational  strain  on  hexagonally  packed  nanotube  array. 
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FIGURE  8.  Calculation  of  Poisson’s  ratio. 


The  C22  stiffness  coefficient  is  obtained  by  applying  the  same  deformation  mode  defined  in  Figure  8, 
but  no  energy  minimization  is  performed  to  induce  a  transverse  Poisson  deformation. 
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5.  Potential  Energy  Calculations.  The  intertube  forces  are  typically  modeled  by  the  Lennard-Jones 
potential  to  represent  van  der  Waals  interactions.  The  Lennard-Jones  or  ‘6-12’  potential  energy  function 
is  given  by 


(5.1) 


v  > 

12 

(  > 

6  “ 

a 

a 

r;i 

|_v  1J ) 

l  1J  ) 

where  q  is  the  depth  of  the  energy  well,  a  is  the  van  der  Waals  radius,  and  r,j  is  the  separation  distance 
between  the  ith  and  jth  atoms  in  a  pair.  (Traditionally,  the  parameters  a  and  'C,  are  referred  to  as  a  and  £, 
which  are  used  in  this  work  for  stress  and  strain  respectively.)  The  r-76  term  represents  the  attractive 

contribution  to  the  van  der  Waals  forces  between  neutral  molecules.  It  implicitly  includes,  via 
parametrization,  permanent  dipole-dipole  interactions,  the  induction  effect  of  permanent  dipoles,  and 
instantaneous  dipole  induced  dipole  interactions  which  are  sometimes  referred  to  as  the  London 
dispersion  forces.  The  other  component  of  the  van  der  Waals  interactions  mimics  the  repulsion  between 
overlapping  electron  clouds  and  is  modeled  by  the  rT 12  term  which  is  short  ranged[18,19]. 


6.  Results  and  Discussion.  In  the  present  study,  a  crystal  of  single-walled  nanotubes  is  analyzed. 
Each  nanotube  has  a  diameter  of  0.94  nm  and  an  achiral  zig-zag  conformation  with  12  graphene  units 
around  the  circumference.  The  carbon-carbon  bond  length  is  prescribed  as  1.42  A  .  Nanotubes  of  this  size 
may  be  considered  rigid  in  the  transverse  direction[12,13].  Therefore,  the  only  degree  of  freedom 
included  in  the  deformation  kinematics  is  the  relative  motion  of  the  nanotube  center,  and  the  only 
contribution  to  the  potential  energy  changes  with  imposed  motion  is  computed  using  the  Lennard-Jones 
potential.  The  parameters  used  in  the  Lennard-Jones  potential  are  q  =  34.0  (K)  and  a  =  0.3406  nm  [13,20]. 
A  comparison  between  predicted  elastic  moduli  using  the  current  direct  method  and  results  obtained  using 
alternate  approaches  is  presented  in  Table  1.  References  13  and  14  utilize  a  lattice  dynamics  approach 
while  Reference  12  is  based  on  a  molecular  dynamic  simulation.  From  the  limited  published  results  it  is 
clear  that  there  is  a  wide  variation  in  predicted  elastic  moduli  for  nanotube  bundles. 


TABLE  1.  Comparison  of  predicted  elastic  constants. 


Elastic  Constant 

Direct  Method 

Reference  13 

Reference  14 

Reference  12 

Bulk  Modulus  K23  (GPa) 

45.8 

42.0 

18.0 

33.6 

Shear  Modulus  G23  (GPa) 

22.5 

5.3* 

— 

— 

Young’s  Modulus  E22  (GPa) 

60.3 

17.0 

— 

— 

Normal  Stiffness  C22  (GPa) 

68.3 

42.0 

78.0 

— 

Poisson  Ratio  Bt, 

0.34 

0.75 

- 

- 

Value  derived  using  relationship  in  Equation  (6). 

Results  generated  using  (7,7)  helical  nanotubes  with  diameter  =  0.94  nm. 
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All  the  results  listed  in  Table  1  for  the  present  analysis  were  computed  independently,  i.e.  none  were 
derived  from  a  subset  of  other  values.  For  a  3-D  solid  exhibiting  hexagonal  symmetry,  the  expression  for 
the  transverse  Young’s  modulus,  E22,  is  given  by[  1 3] 


(6.1) 


(C33  -C32)[(C33  +C32)Cn  —  2C31  ] 

(cJ3c„-c?3) 


Because  the  axial  modulus  of  the  nanotubes  in  the  crystal  is  generally  two  orders  of  magnitude  greater 
than  the  transverse  moduli,  we  may  take  Cn  » C33,C32,C31 .  Using  the  relationship  C23=t>23C22  and 
equating  Poisson  ratios  under  transverse  plane  isotropy  as  u23  =  u32 ,  Equation  (6.1)  reduces  to  a  simple 
relationship  between  C22,  E22  and  n 23 .  Substituting  the  calculated  values  for  the  Young’s  modulus  E22  and 


Poisson’s  ratio B23  into  the  resulting  expression,  the  transverse  normal  stiffness  is  obtained  as 

(6.2)  C,2  =  =  603  2  =  68.2  GPa 

1  — 1)23  1  -  0.34 

which  closely  agrees  with  the  independently  calculated  value  of  C22.  Next,  computing  C23  as  t>23  C22  = 
23.2  GPa  and  applying  the  relations  given  by  Equations  (3.2),  it  is  found  that  the  transverse  plane 
stiffnesses  are  recoverable  from  the  computed  shear  and  bulk  moduli  as 

(6.3)  C22  =  K23  +  G23  =  45.8  +  22.5  =  68.3  GPa 

C23  =  K23  -  G23  -  45.8  -  22.5  =  23.3  GPa 

Thus,  the  computed  elastic  constants  are  completely  self-consistent  for  a  transversely  isotropic  material  as 
required.  The  values  reported  in  Reference  13  (Table  1)  are  also  self-consistent  for  transverse  isotropy 
and  compare  favorably  with  the  present  analysis  for  the  prediction  of  the  bulk  modulus,  K23.  The  results 
from  a  molecular  dynamics  simulation  presented  in  Reference  12  (Table  1)  is  slightly  lower.  The 
calculation  of  the  normal  stiffness,  C22,  which  is  obtained  from  applying  a  similar  deformation  as  that 
used  to  compute  the  bulk  modulus  (applying  only  e22  instead  of  both  e22  and  e33  with  £22  =  e33),  is 
intermediate  between  the  lattice  dynamics  studies  presented  in  Reference  13  and  Reference  14  (Table  1). 
In  Reference  14  nanotubes  of  the  same  radius  were  used  but  with  different  helicity.  However,  the  effect  of 
helicity  is  discussed  in  Reference  13  in  which  it  is  shown  for  several  cases  that  the  configuration  of 
graphene  units  on  the  nanotube  surface  has  a  negligible  effect  for  small  radius  tubes  (  R  <  16  A  ). 

The  C23  value  derived  here  using  G23  and  K23  in  Equation  (13)  or  E22  and  n23  in  Equation  (12)  are  self- 
consistent  and  yield  a  value  of  23  GPa.  The  only  other  available  comparison  for  C23  is  in  Reference  13 
which  presents  a  value  of  32  GPa.  Another  value  not  reported  in  the  literature  is  the  transverse  shear 
modulus  G23  for  which  the  present  analysis  yields  22.5  GPa.  This  value  varies  from  20.2-24.7  GPa  with 
10  %  variation  in  the  Lennard- Jones  a  and  £  parameters.  If  one  applies  the  first  of  the  relations  given  in 
Equations  (3.2)  to  the  data  in  Reference  13,  one  obtains  a  value  for  G23  of  5.3  GPa.  This  shear  modulus  is 
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comparable  to  the  shear  modulus  associated  with  parallel  planes  in  graphite  which  is  experimentally 
measured  as  4.0  GPa[21,22], 

In  the  direct  summation  and  lattice  dynamic  methods,  the  physics  of  cohesion  are  identically 
represented  by  the  same  parameterization  of  the  Lennard-Jones  potential.  Possible  differences  in 
predictions  using  lattice  dynamics  may  be  due  to  the  inherent  integral  averaging  of  force  constants  used  in 
the  lattice  dynamical  matrix  and  the  a  priori  selection  of  interacting  nearest-neighbor  atoms  used  in 
defining  the  primitive  lattice  cell,  both  of  which  are  avoided  in  the  direct  method. 

A  potential  source  of  inaccuracy  affecting  all  methods  is  the  form  of  the  Lennard-Jones  potential 
function  itself.  The  Lennard-Jones  potential  was  originally  developed  for  noble  gases  and  is  known  to 
produce  poor  results  in  other  applications  including  graphite.  It  gives  good  results  for  the  C33  modulus 
(interplanar  separation),  but  the  C44  parallel  plane  shear  modulus  is  under  predicted  by  an  order  of 
magnitude [22].  Alternative  potentials  have  been  proposed[23]  that  yield  accurate  predictions  for  both  the 
transverse  normal  and  shear  moduli  in  graphite.  Additional  study  is  warranted  to  assess  the  spatial 
interactions  of  delocalized  bonds  in  carbon  nanotubes  that  may  be  underestimated  using  a  spherical  atom 
Lennard-Jones  model. 

7.  Conclusions.  A  consistent  method  has  been  formulated  for  and  applied  to  computing  effective 
transverse  mechanical  properties  of  nanotube  crystals.  The  method  is  based  on  specifying  a  unit  cell 
configuration  with  periodic  boundary  conditions,  applying  a  deformation  field  associated  with  a  particular 
strain  mode,  and  utilizing  a  direct  summation  procedure  to  compute  changes  in  potential  energy  from 
which  an  effective  elastic  modulus  may  be  obtained.  The  disparity  between  reported  predictions  of 
mechanical  properties  that  depend  exclusively  on  van  der  Waals  cohesion  and  the  paucity  of  available 
experimental  data  suggest  that  much  additional  investigation  is  warranted  in  this  area.  The  development 
of  a  more  realistic  representation  of  van  der  Waals  interactions  between  nanotube  surfaces  may  be 
required  to  correlate  analytical  predictions  with  future  experimental  measurements  of  nanotube  crystal 
properties. 
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